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We study the prospects for detecting the annihilation products of Dark Matter [DM] in the 
framework of the two highest-resolution numerical simulations currently available, i.e. Via Lactea 
II and Aquarius. We propose a strategy to determine the shape and size of the region around the 
Galactic center that maximizes the probability of observing a DM signal, and we show that although 
the predicted flux can differ by a factor of 10 for a given DM candidate in the two simulation 
setups, the search strategy remains actually unchanged, since it relies on the angular profile of 
the annihilation flux, not on its normalization. We present mock gamma-ray maps that keep into 
account the diffuse emission produced by unresolved halos in the Galaxy, and we estimate that in 
an optimistic DM scenario a few individual clumps can be resolved above the background with the 
Fermi-LAT. Finally we calculate the energy-dependent boost factors for positrons and antiprotons, 
and show that they are always of 0(1), and therefore they cannot lead to the large enhancements 
of the antimatter fluxes required to explain the recent PAMELA, ATIC, Fermi and HESS data. 



I. INTRODUCTION 

Despite many observational and theoretical efforts, the 
nature of one of the main components of the Universe, 
Dark Matter, is still unknown. Indirect evidence based 
on the gravitational interaction between DM and visi- 
ble baryons can constrain the spatial distribution of DM 
particles, assess their cosmic abundance, probe their low 
temperature and non-dissipative nature. However, they 
tell us nothing about the intrinsic properties of DM par- 
ticles. 

This has motivated the search for signals arising from 
the (weak) coupling of the dark sector to ordinary mat- 
ter and radiation, one of the most promising being self- 
annihilation. "Indirect" DM searches are based on the 
search for secondary particles (neutrinos, energetic elec- 
trons, antimatter and 7-ray photons), produced by the 
annihilation or decay of DM particles. If found, these sec- 
ondary particles would constrain the nature of the DM 
particles and provide useful information on its distribu- 
tion (possibly also on small scales). 

Recently a spectacular increase in the positron ratio 
above 10 GeV has been measured by the PAMELA satel- 
lite [1 , while the ATIC balloon-borne experiment has de- 
tected a peak in the total electron+positron flux above 
- 500 GeV [2j. The PAMELA excess has been inter- 
preted in terms of standard astrophysical sources, see e.g. 
Refs. [3] IU [6j [7] , whereas its interpretation as DM anni- 
hilation signal requires unconventional DM particle mod- 
els, see e.g. [BUS], and it is rather severely constrained by 
the absence of an associated flux of IC photons, antipro- 
tons and 7-rays [10] HU Q2J EES HU HI]. Furthermore, 
all DM models with a high annihilation cross section, as 
needed to reproduce the PAMELA data, would heat and 
ionize the baryons in the early universe; the constraints 
that can be set on these models from CMB data do not 
rely on uncertain assumptions on the DM distribution 



in virialized structures, and can therefore be regarded 
as robust and model-independent HO] HZ1 EH]- On the 
other hand the ATIC peak appears to be inconsistent 
with the recent data from the FERMI satellite [IS] and 
the H.E.S.S. Atmospheric Cerenkov Telescope [20] , 

In light of these recent measurements, we discuss here 
the antimatter flux arising from DM annihilations along 
with the associated 7-ray flux. In fact, a 7-ray signal 
from DM annihilation would provide the 'cleanest' evi- 
dence for DM, since photons do not suffer deflection and 
energy losses in the local universe. Besides peculiar spec- 
tral features such as annihilation lines and final state ra- 
diation [2D (55] H31 1211 EH] , an interesting smoking-gun 
for DM would be the detection of many 7-ray sources 
with identical spectra and no counterpart at other wave- 
lengths [10123 [13 [13 [301 [SEISIN 

In addition, a characteristic DM signature may also be 
found in the angular power spectrum of the diffuse 7-ray 
background [39] EE3 EEH HI ESI El S3 US [47]. 

The expected annihilation signal can be factorized into 
a term that only depends on the intrinsic properties of 
DM particles and a term that depends on their spatial 
distribution. For the particle physics properties, we con- 
sider here four benchmark models, representative of the 
most commonly discussed DM candidates, and of the 
models that been invoked to explain the cosmic leptons 
data discussed above. 

The problem of the spatial distribution of DM can be 
tackled in different ways. Analytic methods based on 
the excursion set theory [15] provides a useful, though 
approximate, insight on the evolution of DM halos [4"S] . 
N-body simulations are the best way to study the highly 
nonlinear processes involved in the evolution of substruc- 
tures. Unfortunately, they can only probe a limited range 
of halo masses and scales. The latest numerical simula- 
tions of a Milky Way - size DM halos [50 ] 1ST ] 152 ] [53], 
are able to resolve ~ 100, 000 substructures down to 
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~ 10 Mq at the present epoch. The evolution of micro- 
halos with size close to the free streaming mass can only 
be studied by simulating a small regions at very high 
redshifts [53] . 

As a consequence, modeling the properties of the galac- 
tic subhalos requires aggressive extrapolations which are 
usually performed by means of analytic, Monte Carlo or 
hybrid techniques and therefore are potentially affected 
by large theoretical uncertainties. In this work we rely 
on the results of the Aquarius [52] 153] and Via Lactea 
II |51j numerical simulations. For 7-rays, we then ap- 
ply the hybrid method of [32] to compute the expected 
annihilation flux of 7-ray photons produced within our 
Galaxy. For antimatter, we use the method developed in 
[55] to obtain the boost factor to cosmic-ray fluxes due 
to the presence of subhalos. For the photon flux we also 
need to account for the extragalactic contribution and 
the diffuse Galactic foreground. For the former, we have 
adapted the formalism of [56 sclf-consistently within our 
model. The latter has been modeled through a rescaling 
of the EGRET maps, reducing the flux by 40 % in order 
to account for all new blazars that Fermi-LAT is expected 
to resolve in 1 year. The Galactic foreground has been 
scaled down by 20 %, in order to take into account the 
latest Fermi-LAT data, that do not confirm the presence 
of the so-called "GeV excess" detected by EGRET. The 
antimatter flux is obtained by solving the transport equa- 
tion for high energy positrons and antiprotons produced 
by DM annihilation, and ignoring contributions from as- 
trophysical sources. In this case, the background pro- 
duced by spallation processes is taken from [57] and [5H] . 

The main aim of this work is to assess the reliability 
of this approach, i.e. that modeling the expected DM 
annihilation flux extrapolating the results of state-of-the- 
art numerical simulation provides robust predictions that 
can be used to assess the possibility of detecting the an- 
nihilation signals, of both photons and antimatter, with 
current detectors. 

The paper is organized as follows: In Section [IT] we 
describe our models for DM halos and their substruc- 
tures which contribute to the cosmological part of the 
DM annihilation signal. In Section |III| we introduce our 
particle physics benchmark models that determine am- 
plitude, shape and features of the annihilation spectrum. 
Our model predictions for the 7-ray and antimatter fluxes 
are described in Sections |IV| and \V\ respectively. Finally, 
we discuss our results and conclude in Section IVII 

In this work we adopt the WMAP-5yr [59] flat, ACDM 
model (fl m = 0.26, fl A = 0.74, cr 8 =0.79, n s =0.96, H = 
72 km/s/Mpc). 



II. MODELING DM HALOS AND SUBHALOS 

The recent N-body experiments Aquarius [521 153] and 
Via Lactea II [51] have simulated the DM halo of a Milky- 
Way (MW) -like galaxy in a flat, ACDM model with cos- 
mological parameters consistent, within the errors, with 



those that best-fit the WMAP-lyr and WMAP-5yr data, 
respectively (Aquarius used Q m = 0.25, f^A = 0.75, 
(78=0.9, n s =l, H = 73km/s/Mpc, while Via Lactea 
II used fl m = 0.24, fl A = 0.76, ct 8 =0.88, n s =0.97, 
Ho = 74km/s/Mpc). Thanks to the unprecedented high 
resolution, these simulations were able to resolve sub- 
structures down to masses as small as ~ 10 4 5 M Q ( 
Aquarius) and ~ 10 5 M Q (Via Lactea II), to character- 
ize their inner structure, to trace their spatial distribu- 
tion within the main halo and to model the dependence 
of their shape parameter (the concentration) on the dis- 
tance from the Galactic Center (GC). For the Via Lactea 
II we will consider subhalos selected by mass as in [BD] 
rather than by peak circular velocity as in [51]. We have 
checked that using either subhalo parameterizations does 
not significantly affect our predictions for the detectabil- 
ity of the annihilation signal. 

Following the results of numerical simulations, the DM 
distribution in the MW halo consists of two separate 
phases: a smoothly distributed component (main halo) 
and a clumpy component made of virialized substructures 
(subhalos). We therefore ignore the presence of caustics, 
streams and all other possible inhomogeneities that do 
not correspond to virialized structures. 

The smooth density profile of the MW DM halo can be 
modeled by a Navarro, Frenk and White (NFW) profile 
in the case of Via Lactea II : 
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and by a shallower Einasto profile with a 
Aquarius: 
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where R is the distance from the GC. The best fitting 
values of the scale density, p si and the scale radius, r s , 
are listed in Table [I] The density profile of the smooth 
components are shown in Fig. [I] 

Being Mh and M sub the main halo and subhalo masses, 
the joint spatial and mass subhalo distribution is best 
fitted by the following relations [5T1 152"1 1551 |S0] : 
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the case of Via Lactea II and 
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with a — 0.678 for Aquarius. Units are M^kpc — . The 
scale radii of the subhalo distributions R a are listed in 
Table [I] The normalizations A s f l (Mf l ), shown in Table [T] 
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FIG. 1: Density profiles of the smooth halo of the MW for 
the Via Lactea II and Aquarius cases. Their normalization is 
described in the text. 



are also chosen to match the results of N-body simula- 
tions. In the case of Via Lactea II, we impose that 10 % of 
the MW mass, Mh, is in virialized structures with masses 
in the range [10 -5 Mh, 10 -2 Mh}. In the case of Aquarius, 
we require that 13.2 % of Mh is concentrated in subhalos 
with masses in the range [1.8 x lO~ 8 M h , 10~ 2 M h ]. 

We note that the MW mass in both simulations agrees, 
within the errors, with the recent observational estimates 
of [61] based on the so-called Timing Argument [62] . 

We also note that in the two numerical experiments the 
smooth halo mass density profile and the subhalo number 
density profile have the same overall shape (NFW in Via 
Lactea //and Einasto in Aquarius ), but are characterized 
by different shape parameters. The total number density 
of substructures with masses above 10 6 M© in Via Lactea 
II is ~ 20 times larger than in Aquarius. However, the 
dumpiness fraction /, defined as the ratio of the mass 
virialized in substructures to the total virial mass of the 
halo, turns out to be similar. These differences arise from 
the different background cosmology assumed in the two 
simulations and from the way in which subhalos were 
selected. 

A word of caution is required for the subhalo distribu- 
tion near the GC (e.g. [63] ). Current numerical exper- 
iments are able to resolve only a few subhalos within 8 
kpc from the GC so that some hypothesis is required to 
model their distribution in the central region of the MW. 
In this work we make the educated guess by extrapolating 
at small separation the subhalo number density profile 
measured at larger radii. In fact, tidal effects may dis- 
rupt subhalos in the central regions of the Galaxy, which 
severely depletes the subhalo population. To account for 
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FIG. 2: Maximum subhalo mass that can be found at distance 
R from the GC, according to the Roche criterion used in this 
paper. 



this effect we adopt the Roche criterion [64 : a subhalo 
is destroyed when its scale radius r s is larger than the 
tidal radius, i.e. the radius at which the tidal forces of 
the host potential equal the self-gravity of the subhalo: 

-(Mi?)" 3 * (5) 

where R is the distance from the subhalo center, M su b 
the subhalo mass and Mh the host halo mass. Clearly, 
the amplitude of the effect depends on the subhalo 
mass and its distance from the GC. In Fig. [2] we plot 
the largest mass A/™g x of a subhalo that survive tidal 
disruption as a function of R. Subhalos are shown to 
be almost completely disrupted within R ~ 2 kpc. As 
shown by [57], the effect on the 7-ray flux is negligible. 
This turns out to be the case also for the antimatter flux 
as we will show, and as already pointed out by [65]. We 
will discuss this point in greater detail in Section [VC 1| 

A very informative quantity to describe subhalos is 
their concentration parameter, defined as the ratio be- 
tween 7' 2 oo (the radius which encloses an average density 
equal to 200 times the critical density of the universe) and 
the scale radius: C200 = r 200 l r s- In general this quantity 
is not constant but depends on the subhalo mass and on 
the distance from the GC. Following the numerical results 
of [5TJ [55] [53] |5D] we can parameterize these dependences 
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Via Lactea 


77 Aquarius 


Rvir [kpc] 


402 


433 


M h [M ] 


1.93 • 10 12 


2.5 • 10 12 


r s [kpc] 


21 


20 


p s [1O 6 M kpc- 3 ] 


3.7 


2.3 


A 8h [Mg 1 kpc" 3 ] 
fl a [kpc] 


1.7 x 10 4 
21 


25.86 
199 


N cl 


2.79 • 10 ie 


1.17- 10 15 


M|'(< 7? mr )[M ] 


1.05 ■ 10 12 


4.33 ■ 10 11 


/c? (< Rvir) 


0.53 


0.17 



TABLE I: Characteristics values for the smooth and clumpy 
components of the DM distribution modeled after the Via 
Lactea II and Aquarius results. Rvir = virial radius, i.e. 
the radius within which the numerical simulations define the 
halo. Mh = MW mass. r s = scale radius of the smooth DM 
distribution. p s = scale density of the smooth DM distribu- 
tion. A s h = normalization of the subhalo distribution. R a = 
scale radius of the subhalo distribution. N c i = total number 
of subhalos. A7*°' = mass in subclumps. /^ ot = dumpiness 
fraction. 
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Via Lactea II Aquarius 


CtR 


0.286 


0.237 


CI 


119.75 


232.15 


C2 


-85.16 


-181.74 


aci 


0.012 


0.0146 


ac2 


0.0026 


0.008 



TABLE II: Parameters for the fit to the concentration param- 
eter of subhalos. 



as follows: 



c 2 oo{M sab ,R) = ( - — J x 

Ht?r+-(t?r 

The best fitting parameters for the Via Lactea II and 
Aquarius simulations are listed in in Table [TTJ In Fig. [3] 
we have plotted the mass dependence of the concentra- 
tion parameter at the virial radius (7? — R V i r ) that can be 
thought of as the concentration parameters of subhalos 
located at the edge of the simulated MW-like halo, i.e. 
of fields halos. We have also plotted the concentration 
parameters computed at the Sun position (7? = 8 kpc), 
which provides additional information on the potential 
antimatter yield by featuring the local subhalo prop- 
erties. From the plot we notice that subhalos in the 
Aquarius experiment are more concentrated than in Via 
Lactea II at all masses; a discrepancy that reflects the 
larger power spectrum normalization (er 8 ) assumed in the 
Aquarius experiment. 



FIG. 3: Concentration parameter as a function of the halo 
mass for field halos, as found in Via Lactea II and Aquanus, 
computed at a the virial radius Rvir and at the Earth-GC 
distance R = 8 kpc. 



III. PARTICLE PHYSICS BENCHMARKS 

In order to study the dependence of the results on par- 
ticle physics parameters, we show the results relative to 
four different benchmark models. 

Benchmark A is representative of a class of mod- 
els in the supersymmetric (SUSY) parameter space, that 
annihilate predominantly to bb. In order to maximize 
the annihilation flux, we chose a light neutralino mass 
{Q)n x = 40 GeV. 

Benchmark B is also representative of a class of 
2 \SUSY models. The DM particle mass is in this case 
Jm x — 100 GeV, thus allowing annihilation to W + W~ , 
that is assumed to constitute the dominant annihilation 
channel. 

Benchmark C provides a "minimal" solution to the 
rising positron ratio measured by PAMELA, without at- 
tempting to address higher energy (ATIC and Fermi) 
data. The mass is in this case m x = 100 GeV, thus 
barely above the PAMELA energy range, and the lead- 
ing annihilation channel e + e~. 

Benchmark D, finally, represents a class of candi- 
dates that attempt to explain at the same time the cos- 
mic lcpton data up to TeV energies. We have adopted, 
as e.g. in Ref.2008PhRvD..78j3520B, m x = 2000 GeV, 
and annihilation to t + t _ . 

We have adopted in all cases a thermal annihilation 
cross section (av) = 3 x 10 _26 cm 3 s _1 . The parameters 
of the four benchmark models are summarized in Table 

urn 
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model m x [GeV] final state 



A 


40 


bb 


B 


100 


w+w~ 


C 


100 


e + e~ 


D 


2000 


r+r- 



TABLE III: Benchmark particle physics models. The annihi- 
lation cross-section is (av) — 3 x 10 -26 cm 3 s _1 . 



IV. GAMMA-RAYS 



The expected 7-ray flux from DM annihilation, $ 7 , can 
be factorized in two terms that depend on the properties 
of DM particle, Qpp, and on their spatial distribution 
along the line-of-sight, <&i os : 



<$> 1 (m x ,E 1 ,M,r,d,e) = $ P p(m x ,E^) x ® los (M,r,d,6) 

(7) 

where m x is the DM particle mass, M the DM halo mass, 
r the position inside the halo, d the distance from the 
observer and 9 the angular resolution of the instrument 
(in the case of Fermi, for energies above ~ 1 GeV, one 
has 6 ~ 0.1°). 

The term Qpp describes the number of photons yielded 
in a single annihilation, and can be written as: 



fbp P [M,E~) 



An 2ml J Ei y dE i h 



(8) 



where / is the final state, Bf is the branching ratio, and 
(av) denotes the thermal annihilation cross section which 
reproduces the observed cosmological abundance today: 
(av) = 3 x 10 _26 cm 3 s -1 . dNj/dE^ is the differential 
annihilation photon spectrum that we take from [66 . 

The term <&; os represents the number of annihilation 
events along the line-of-sight. It is obtained by integrat- 
ing the square of the DM mass density: 



$, os (Af,r,d,0) = 

p 2 DM (M, c,r(A,y>,fl,0)) 
d 2 



d9d(j) / dX 



J(x,y,z\\,9, 



(9) 



where Af2 is the angular resolution of the instrument, J 
is the Jacobian determinant, d = A and c the concentra- 
tion parameter. In the case of the smooth halo of the 
MW, M = M h and c is fixed by the output of N-body 
simulations, while for the subhalos M — M su f, and c is 
a function of mass and position: c — C2oo(M su b, R), as 
defined in Eq. ^ . 

The 7-ray annihilation flux receives contribution from 
three different sources that we model separately: 

• the DM smoothly distributed in the MW halo 

• the DM within Galactic subhalos 

• the DM in extragalactic halos and their substruc- 
ture. 



To compute the first contribution we simply consider 
the mass density profile in Eq. or Eq. ([2| in the 
integral Eq. 

To compute the second contribution one would need 
to consider all ~ 10 18 substructures down to the cutoff 
mass of 10 ~ 6 Mq, which is unfeasible. One possibility 
would be to integrate the subhalo distribution (Eq. [3] or 
Eq. [4| . This represents the mean annihilation flux while 
one of the scopes of this work is to assess the possibil- 
ity of detecting isolated DM subhalos. To circumvent 
the problem we follow the hybrid approach of [32], i.e. 
we first compute the mean flux and then we use Monte 
Carlo techniques to account for the closest and brightest 
subhalos that one may hope to detect as isolated sources. 

More specifically, the mean flux along the line-of-sight 
is obtained by integrating the following expression: 



<C s b (A4,d>,Aft)cx 



dM,. 



dc 



dX[p sh (M h , M sub , R)P(c)<S> los 



(10) 



where $; os represents the contribution from each sub- 
halo computed by integrating Eq. |9| in the range [d — 
Tudi d+r t id\- This quantity is convolved with the subhalo 
distribution (Eq. [3] or Eq.Q. 

P(c) represents the probability distribution function of 
the concentration parameter. Following |67j . we model 
it as a lognormal distribution with dispersion a c = 0.14 
and mean value c: 



P(c,c) = 



1 



ln(c)-ln(c) 



2ira c c 



(11) 



To model the fluctuations over the mean flux we com- 
pute the annihilation flux from the nearest and bright- 
est subhalos in 10 independent Monte Carlo realizations. 
For this purpose we consider only those subhalos whose 
distance from the Sun is the maximum between: 1) the 
radius of the sphere centered on the sun within which lie 
about 500 halos, and 2) the distance at which the photon 
flux from a subhalo drops below the value of the aver- 
age flux at the anticenter. The results may depend on 
the actual number of individual subhalos in the Monte 
Carlo simulations. To check the robustness of the re- 
sults [32 have performed a number of convergence tests 
in which they demonstrated that increasing the number 
of individual halos does not change the estimate of their 
detectability. The scatter among the results is obtained 
from the different Monte Carlo realizations of the subhalo 
population. 

Finally, to compute the extragalactic contribution to 
the annihilation flux we have used the formalism of [56J, 
modified to account for the mass and radial dependence 
of the subhalo concentrations. Here we outline the main 
steps of the computation. A detailed description can be 
found in [T2l. 
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FIG. 4: Left and central panel: Contribution to the 7-ray flux from DM annihilation, as a function of the angle of view from 
the GC, of the MW diffuse (smooth red curve), MW subhalo (irregular blue) and extragalactic (horizontal green) components 
for the Aquarius (left panel) and Via Lactea II (central panel) setup. The right panel shows the sum of all contributions for 
the Aquarius (red curve) and Via Lactea II (green) setup. 



The contribution to the annihilation flux of isolated 
field halos can be obtained by integrating the contribu- 
tion of a single object dMj/dE over the halo mass func- 
tion 1 68 : 



Good fits are 



dn 
dM 



1 ' ; M 2 Jy 'dlosM' y ' 



where p cr is the critical density, f2o,m is the mass density 
parameter, v — S °^} , cr{M) is the rms density fluctua- 
tion on the mass scale M and S sc represents the critical 
density for spherical collapse [55] EHJ [70] • 
In addition, the attenuation for pair-production along 
the line-of-sight is described by introducing an absorp- 
tion coefficient e~ T ( z, - E °) [71], where E is the energy of 
the photon as it is observed today. As a result the ex- 
pression for the extragalactic DM 7-ray flux is given by 



with 



A 2 (z) = J dM 



v(z,M)f(v(z,M)) 
Ma(M) 



da 



dM 



MA 2 M (z,M) 
(14) 

A 2 (z) is a quantity which describes the boost to the 
isotropic 7-ray background due to the existence of viri- 
alized DM halos, g(z) is a redshift dependent quantity 
related to the expansion history of the universe, and 
A M cx $i os describes the single halo contribution. 

To account for the fact that a fraction f(M) l ) of the 
main halo of mass is located on virializcd subhalos we 
have modified the single halo contribution as follows: 

M h A 2 M (M h ) -> (l-f(M h )) 2 M h A 2 M (M h ) + A 2 MiSub (M h ) 

(15) 

The function /(M/J was modeled after the numerical 
simulations and the extrapolation down to 10~ 6 M Q . 



f(M h ) = 1.45 x 10- 2 [log(M h ) + 9.21] 



in the case of Via Lactea II and 



f(M h ) = 0.285Af^ 01 [0.63M^ 1 - 0.25] 
in the case of Aquarius. 



To compute the extragalactic flux we integrate Eq. ( 13 1 



numerically taking into account both the main halo and 
subhalo contributions. The latter is provided by the 
Am, sub term which we also evaluate numerically taking 
into account the subhalo distribution (Eq. [3] or Eq. [I| , 
the probability distribution function of the concentration 



parameters (Eq. Ill and their dependence on mass an 



distance (Eq. 6), as follows: 



A 2 Msub (z,M) cx / dM sub I dR I dc 



(16) 

M sub JR 

Psh (M h , R)P(c(M sub , 2, R))$i os (M sub ) 

The left and central panels of Fig. [4]show the predicted 
contribution of the smooth MW halo, MW subhalos and 
extragalactic halos to the 7-ray flux from DM annihila- 
tion integrated above 3 GeV, as a function of the angle 
of view ip from the GC in the case of Aquarius and Via 
Lactea II, respectively. The contribution from individual 
subhalos is computed by averaging over the 10 Monte 
Carlo realizations. 

In the central part the annihilation signal is dominated 
by the MW smooth component in both Aquarius (at ip < 
50°) and Via Lactea II (at i\> < 20°). Away from the GC 
the dominant contribution is provided by DM subhalos in 
the Via Lactea II case and by extragalactic background in 
the Aquarius case. The different behaviors simply reflect 
the different amount of Galactic substructures in the two 
simulations. 



FIG. 5: Full-sky map, in Galactic coordinates, of the number of photons (above 3 GeV) produced by DM annihilation (bench- 
mark A). The left (right) panel shows the predicted flux in the Aquarius (Via Lactea IT) setup. 




FIG. 6: Same as Fig. [5j but with the two simulation setups 
fraction of mass in substructures. 

We note that our predictions satisfy the observational 
constraint represented by the diffuse Galactic signal. In- 
deed the mean diffuse Galactic flux above 3 GeV mea- 
sured by Fermi (~ 5.3 x 1CP 7 ph cm -2 s _1 sr _1 ) is safely 
above the expected annihilation signal in both simulation 
setups. 

The flux in the innermost regions is higher for the 
Aquarius simulation, also plotted in the full sky maps 
of Fig. [5] where we show the total annihilation flux (MW 
smooth + galactic subhalos + extragalactic halos and 
subhalos). The fact that the annihilation signal at the 
GC is stronger for the Aquarius simulation is due to 
the combined effect of: i) the larger fraction of mass in 
substructures in the Via Lactea II setup, that reduces 
the annihilation flux from the smooth component by a 
factor (I — /^ 9 ) 2 /(f — f^ L2 ) 2 ~ 3; 11) the smaller lo- 
cal density, i.e. DM density in the solar neighborhood, 
in the Via Lactea II setup, that reduces the the anni- 
hilation flux from the smooth component by a factor 
(p©7p© L2 ) 2 = (0.48/0.2) 2 = 5.8 and Hi) the smaller 
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rescaled to the same local density, same total mass and same 



total mass in the MW in the Via Lactea //simulation. 

If one rescales the Via Lactea II and Aquarius setup to 
match the most recent determinations of the local DM 
density p G = 0.385 GeV/cm 3 [H [73], and if the same 
subhalo mass fraction / c ; = 0.18 is adopted for both sim- 
ulations, then the annihilation maps look almost identi- 
cal, as shown in Fig. [6] 



A. Experimental detectability 

In order to assess the detectability of the 7-ray annihi- 
lation flux with the most promising experiment currently 
available, i.e. the Fermi-LAT satellite, we have to specify 
what is our signal and what is our background or noise. 

If we are interested in finding a signal above the as- 
trophysical backgrounds, the signal is contributed by the 
sum of all the aforementioned components of the anni- 
hilation flux (MW smooth mass distribution + galactic 
subhalos + extragalactic halos and subhalos). We fo- 



FIG. 7: Sensitivity map, in galactic coordinates, for the Aquarius and Via Lactea //setups. The signal is as in fig.??, while 
the background is obtained through a suitable rescaling of the EGRET maps (see text for further details). 



cus on photons with energies larger than 3 GeV collected 
over 1 year of on-axis observation, which corresponds to 
about 5 years of data taking with Fermi, and we assume 
an effective detection area of 10 4 cm 2 . The background 
or noise is contributed by the diffuse Galactic foreground 
and the unresolved Extragalactic background. To model 
the Galactic foreground we have scaled down the signal 
measured by EGRET by 20% to take into account the 
fact that Fermi-LAT preliminary result did not confirm 
the EGRET GeV excess. The extragalactic background 
of Fermi can be estimated by rescaling the EGRET back- 
ground down by 40 % to account for the blazars that will 
be resolved by Fermi. The expected sensitivity is sim- 

JV 7 

ply given by a = ^/ N £^ ^ wh ere N s ] gnal is the to- 
tal number of annihilation photons and -^background are 
the photons contributed by the astrophysical background 
and foreground. The left and right panels of Fig. [7] show 
the resulting sensitivity maps in Galactic coordinates for 
the Via Lactea II and the Aquarius simulations, respec- 
tively. The sensitivity maps have a sharp peak near the 
GC, as expected. Iso-sensitivity contours are not sym- 
metric but have a characteristic 8-shape around the GC 
due to the disk-like Astrophysical Galactic Foreground 
that contribute to the noise term. The sensitivity maps 
after the rescaling procedure adopted to produce the two 
maps in Fig(6] would look very similar to the left panel of 
Fig. 

This procedure can be useful to estimate the regions 
that optimize the signal-to-noise in DM searches. In fact, 
once the Fermi data become available, one can just take 
the DM template presented here, and divide the actual 
Fermi data by such template. One can then calculate 
the signal-to-noise ratio, S/N, for the iso-flux contours, 
and determine the size and the shape of the region that 
maximizes the S/N. 

Disentangling the annihilation signal from the astro- 
physical one near the GC might however be difficult due 
to the presence of a strong astrophysical background. An 



alternative strategy is to look for individual subhalos, 
i.e. isolated bright spots in the 7-ray sky. In this case 
the signal is given by the annihilation photons produced 
in the nearest and 7-ray brightest subhalos in our Monte 
Carlo realizations. The noise is contributed by all 
remaining sources of 7 ray photons, including those from 
DM annihilation (MW smooth + extragalactic). The 

sensitivity is in this case given by a = " emI = where 

V -^back+forc 

signal = JV 7ubhalos,MC are the annihilation photons 
produced within the resolved subhalos and N^ ck+foTc — 
/V 7 +/V 7 +/V 7 + /V 7 

MW, smooth subhalos, average extragalactic background 

includes contributions from the smooth MW halo, un- 
resolved subhalos, diffuse extragalactic background and 
Galactic foreground, respectively. 

In Table II VI we list the number of 3 a and 5 a 
detections expected with an exposure of 1 yr on-axis 
with the Fermi-LAT. Both the Aquarius and the Via 
Lactea //cases are considered. The number of detectable 
subhalos is small but significantly different from zero 
in benchmark models A and B. No individual subhalos 
are expected to be detected if the case of models C and D. 

Fig. [8]shows the distribution of the subhalos detectable 
at 3 and 5 a for the benchmark model A, as a func- 
tion of the subhalos mass. The symbols and error bars 
represents the mean and the scatter over the 10 Monte 
Carlo realizations. All detectable subhalos have masses 
above 1O 5 M0, in some cases comparable to the estimated 
masses of the local dwarf galaxies, suggesting that these 
DM-dominated objects may indeed be good targets for 
DM indirect detections. Furthermore, the results show 
that the only detectable halos are those already resolved 
in the Aquarius and Via Lactea //simulations, i.e. the 
results presented here are independent on the aggressive 
extrapolations required to model the properties of low- 
mass subhalos. 
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FIG. 8: Number of halos that can be detected at 3 a with an 
exposure of 1 year on-axis (corresponding to approximately 
5 years of data taking) with the Fermi-LAT, in the case of 
the Via Lactea II and Aquarius simulations and the particle 
physics benchmark model A, as a function of the halo mass. 
Errors represent the scatter over 10 Monte Carlo realizations. 



model 


VLII 3 a 


VLII 5 cr 


Aq 3 a 


Aq 5 a 


A 


9.2 ±2.6 


4.1 ± 1.3 


13.5 ±2.5 


7.3 ± 2.4 


B 


3.1 ± 1.1 


1.4 ±0.8 


6.2 ±2.1 


2.9 ± 1.4 



TABLE IV: Number of subhalos detectable at 3 and 5 a with 
the Fermi-LAT satellite in 5 years of data taking, computed 
for the benchmark models A and B. Errors represent the scat- 
ter over 10 Monte Carlo realizations. 



V. ANTIMATTER 

Antimatter cosmic rays [CRs] have long been consid- 
ered as a potential tracers of DM annihilation because 
they are barely produced in standard astrophysical pro- 
cesses (T4j . Indeed, most of the standard contributions 
are expected to be of secondary origin, i.e. produced by 
nuclear interactions of standard CRs (protons and light 
nuclei) with the interstellar gas (hydrogen and helium). 
This picture is essentially true for antiprotons, because 
their propagation scale is rather large compared to the 
spatial fluctuations of the interstellar medium and CR 
sources, but should be taken with caution for positrons, 
since their high energy component is strongly sensitive 
to time and spatial fluctuations of the local environment. 
Moreover, some astrophysical sources, like pulsars, are 
known to have the capability to produce positrons from 
pair creations in strong magnetic fields (sec an early dis- 
cussion on this in |75j). Because there are such sources 
in abundance within the kpc scale around the Earth, one 



can expect them to be significant contributors to the high 
energy positron budget j4]. Nevertheless, many theoret- 
ical as well as observational uncertainties still affect the 
modeling of most of astrophysical sources. It is conse- 
quently important to scrutinize the potential imprints 
that DM annihilation could provide in the antimatter 
budget, in addition to those in 7-rays. In this section, we 
will shortly review the transport of CRs in the Galaxy, 
and then focus on the positron and antiproton fluxes at 
Earth that DM annihilation may generate. To model the 
sources of the annihilation products we will again use the 
highest resolution N-body simulations discussed in the 
previous sections, including the effect of subhalos. For 
more details on the semi-analytical method used to fol- 
low the antimatter transport within the Galaxy that will 
be used in the following, we refer the reader to [531155] , 



A. High energy positron and antiproton transport 

Antimatter CRs in the GeV-TeV energy range, like 
other charged cosmic rays, diffuse on the inhomogencitics 
of the Galactic magnetic field. Because those inhomo- 
geneities are not fully confined in the Galactic disk, cos- 
mic rays can pervade beyond this tiny region and dif- 
fuse away up to a few kpc upwards and downwards (see 
e.g. [80] for a pedagogical insight). This turbulent vol- 
ume actually defines the diffusion zone inside which CRs 
are confined, and beyond which they escape forever; it 
can be featured like a disk-like slab with radial and ver- 
tical extensions, -R s i a b and L, respectively. In the fol- 
lowing, we will fix i? s i a b = 20 kpc, and L = 4 kpc |77j . 
CRs experience different processes during their journey, 
depending on their nature: in addition to spatial diffu- 
sion, (anti)nuclei will be mostly affected by convection 
and spallation in the interstellar medium localized in the 
disk, processes which are more efficient at low energy, 
while electrons and positrons will have their transport 
dominated by energy losses above a few GeV. Formally, 
given a source Q, all CRs obey the same continuity equa- 
tion [75] . 

VJ = d^ + d E J E -TN= Q(x,E,t) , (17) 

where M = dn/dE is the CR density per unit energy, 
is the space-time current, J E the energy current and 
r stands for a destruction rate (spallation for nuclei). 
The time current is merely the CR density J 1 = N , 
while the spatial current is reminiscent of the Fick law, 
and accounts for the spatial diffusion and convection 

J s = {K X (K)\7 + V c }aT. The spatial diffusion coef- 
ficient, describing the stochastic bouncing interactions 
with the magnetic inhomogencitics, is usually parame- 
terized as K X (1Z) — K (TZ/1 GV) S where the dependence 
on the CR rigidity TZ = p/Z is explicit. The energy 
current carries the energy loss and reacceleration terms 
J E = {dE/dt + K p {E)d E }M. 



10 



Though all CR species obey the same transport equa- 
tion, some of the processes mentioned above will be neg- 
ligible in the GeV-TeV energy range, depending on the 
species. Beside spatial diffusion, antiprotons will mostly 
suffer spallations and convection, but almost not energy 
losses. On the contrary, energy losses will dominate the 
positron transport, mainly due to inverse Compton scat- 
tering with the CMB photons and the interstellar radi- 
ation fields, and to synchrotron losses with the Galactic 
magnetic field. In both cases, reacceleration is negligi- 
ble above a few GeV [55J E7]- Neglecting the irrele- 
vant terms and assuming steady state (d t N — 0), the 
Green functions — or propagators — associated with 
the transport equation can be derived analytically for 
antiprotons and positrons. We refer the reader to |65j 
for the detailed expressions of the propagators. In the 
following, Q will denote the Green function, such that 
VG(x,E <- x s ,E s ) = S 3 (x-x s )8(E-E s ); S indexes the 
source quantity. If one forgets about the spatial bound- 
ary conditions, it is useful to remind that the antipro- 
ton propagator is oc S(E — Es)/(4:irK('R)r), while the 
positron propagator is cx 1/ X 3 exp(—(r / X(E)) 2 ) , where r 
is the source distance and X(E) is the energy-dependent 
propagation scale. More details on propagation can be 
found in e.g. [75J E3 EQl EH [82] . 

Most of the propagation parameters can be constrained 
from measurements of the ratios of secondary to par- 
ent primary species, and we will use the median set 
derived in [53]: K = 1.12 • 10- 2 kpc 2 Myr- 1 , 8 = 0.7 
and V c — 12 km/ 's. The energy loss rate ascribed to 
positrons will be dE/dt = —b E 2 , with fe — ( r ioss x 
lGeV) -1 = 10 _16 s _1 GeV _1 , which is a reasonable ap- 
proximation accounting for the inverse Compton loss 
on CMB, starlight and dust radiation and for the syn- 
chrotron loss 58, 80, 84 . However, we remind that there 
are theoretical uncertainties and degeneracies among 
those propagation parameters [77], so that predictions 
for primaries and secondaries may vary by large factors 
(see |83j for primary antiprotons and [65 , 85 for primary 
positrons) . 

An important consequence of the differences in the 
propagation histories among species is that the corre- 
sponding characteristic propagation scales also differ. 
For antiprotons, the propagation scale is set by the spa- 
tial current, X p = K(K)/V C ~ 1 (T/GeV) <5 kpc, and in- 
creases with energy. On the contrary, positrons have 
their propagation scale set by energy losses, A e + cx 

{-/ dEK(E)/(dE/dt)} 1/2 ~ 7(£/GeV) 5 - 1 kpc, which 
strongly decreases with energy (see e.g. Fig. 2 of [86 J. 
Therefore, low/high energy antiprotons/positrons ob- 
served at Earth must originate from the very local en- 
vironment. For instance, a ~100 GeV positron has a 
probability to originate from regions farther than ^2 kpc 
Gaussianly suppressed. 

The CR flux at the Earth for any source Q(xs,Es) 



reads therefore 

4>{E) = f d 3 x s x 

J slab 

p CO 

/ dE s g(x Q ,E «- 2 S , E s ) Q(2 S , E s ) ■ 

JE 

(18) 

It proves convenient to separate the spatial and en- 
ergy parts by writing the source term as Q(xs,Es) = 
Qo(2s)j§^(Es), where Qo carries the units of inverse 
volume and inverse time, and 4£f- is merely the num- 
ber of CRs injected within dEs- The expression above 
thereby simplifies as 

<f>(E) = f d 3 x s G(2 e «- 2 S , E) Qo(xs) , (19) 

J slab 

where we have defined 

~ f°° dN 

G(x Q «- 2 S , E) = / dE s g(x Q ,E «- x s ,E s ) — (20) 
Je dE s 

We remind that for antiprotons, G °c 8(E — Es) since 
we neglect energy losses and reacceleration, so that the 
injected CR spectrum and the propagator can be fac- 
torized. This is not true for positrons, for which energy 
losses are important and prevent from such a simplifica- 
tion. 

Beside the primary signals due to DM annihilation, 
that we will discuss below, one should also be aware of the 
backgrounds and their theoretical uncertainties. Since we 
disregard here the conventional astrophysical sources of 
CRs, our positron and antiproton backgrounds are those 
secondaries resulting from spallation processes of cosmic 
protons and nuclei with the interstellar matter located in 
the disk. We refer the reader to [57] and [58] for thorough 
discussions on those secondary components and related 
theoretical uncertainties. Throughout the paper, we will 
use the median secondary backgrounds derived in those 
references. 



B. Smooth and clumpy DM contributions: boost 
factors 

The fact that the DM spatial distribution is not 
smooth, but actually fluctuates due to the presence of 
subhalos, leads to local fluctuations in the annihilation 
rate [87]. Formally, any flux estimated from a smoothly 
averaged DM profile should therefore be enhanced by 
a factor (p 2 lm )y cr /(/ dm) 2 / to account for those fluctu- 
ations, the average being performed in a volume V CI that 
depends on the CR propagation scale. Such an enhance- 
ment must be, therefore, quite different from what has 
been previously discussed for 7-rays, simply because the 
averaging volume for the latter is the resolution cone car- 
ried by the line of sight instead of a propagation volume. 

Using the notations defined in Section |V A| the anti- 
matter flux at the Earth originating from the annihilation 
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of a single, smoothly distributed DM component is given 
by the following expression 



(3c 
Air 



S 



d 3 xsG(x e ^xs,E) 

lab V P& 



(21) 



where Q already contains the injected spectrum dN/dEs, 
S = ((av) /2)(p Q /m x ) 2 and p Q = p sm (x Q ) is the DM 
density near the Sun inferred from the smooth profile 
p sm in Eqs. and Q. 

Besides, we need to determine the average contribution 
due to the population of subhalos. The number density 
of subhalos of mass M sub and at position R given by 
Eqs. ^ and @ can be used to obtain the corresponding 
normalized probability function as follows: 



p sh (M sub ,R) = N t 



lot dV v (R)dV m (M S ub,R) 
dV dM sub 



(22) 



where N*°^ is the total number of subhalos. This ex- 
pression is such that the product of dVs corresponds to 
a normalized probability function 



dVt 



mw dV 



dm 



dV m (m,R) 



dm 



(23) 



Defining such a probability function will allow us to treat 
each quantity related to a single subhalo as a stochastic 
variable [53]. Notice that though the upper bound in 
the integral of the mass distribution is fixed to m ma x, 
the mass distribution itself is in fact a function of R. 
Such a dependence arises because of the tidal disruption 
of subhalos that is implemented in the present analysis 
according to the Roche criterion (Eq. [5]) that we have 
discussed in Section ITT1 

Before inferring the overall subhalo flux, we need to 
define the luminosity of a single object. Given a subhalo 
of inner profile p su b DM , the corresponding annihilation 
rate will be proportional to the annihilation volume 



Zsub(M sub ,R)= 4n / drr 2 

JVgub 



Psub 



DM (r,M sub ,R) 



(24) 



which is the volume that would be necessary to obtain 
the whole subhalo luminosity from the local DM density 
Po associated with the DM setup of interest. As de- 
tailed in previous sections, the concentration parameter 
c(M sub , R), which characterizes the shape of the inner 
density profile, depends on the subhalo mass M sub and 
position R in the Galaxy. 



Therefore, the CR flux for a subhalo population reads 

&E) = ^SN^ f d^M^E)^ r^ dm Ub( m ,R) dV ^ R) 

47r ./slab dV J mm . n 

^SN& [ d ^ s g {Ss ,E)^P-(U b ) Ma jR) 

J slab aV 



dm 



4tt 

= ^SN™ (Q(x s , E) (Uh) Msub (R)) slah = (0 sub ) . (25) 

( x )m/v me a n s an average of the variable x according to the mass/spatial distribution, respectively. By writing this 
equation, we made, as mentioned above, the implicit assumption that a subhalo can be treated as a stochastic point- 
like source. This assumption is valid while the typical propagation scale is larger than the scale radius of any subhalo, 
which is fully the case here: in the Aquarius setup, a 1O~ 6 /1O 6 M0 subhalo has a scale radius of ~ 10~ 6 /3 • 1CP 2 kpc, 
respectively. We obviously recover the same result as in [5S] that the total subhalo flux is given by the total number 
of subhalos times the mean flux from a single subhalo, consistently with the stochastic treatment. Nonetheless, the 
spatial average cannot be separated from the mass average anymore in the present study. This comes not only from 
the full implementation of tidal effects, but also from the spatial dependence of the concentration parameter (the 
two effects are related). Consequently, the average luminosity of a subhalo, (£sub)jv/ 6 > does explicitly depend on 
the location in the Galaxy, even when the spatial disruption is neglected. We illustrate the radial dependence of 
(£sub)j\,f b (R) in Fig. [9j where we explicitly show the effect of the tidal disruption: the average luminosity of a subhalo 
is quickly turned off m the central part of the Galaxy when Roche criterion is applied. 

Because the flux derived above is in fact a mean quantity, it can be associated with a statistical variance. The 
relative fluctuation of the antimatter CR flux is given by 



'M gub y "slab 



This relative fluctuation, associated with the whole pop- 



ulation of subhalos, is linked to the fluctuation of the 
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Galactocentric radius R [kpc] 



FIG. 9: Mean value of the subhalo annihilation volume de- 



fined in Eq. (241 as a function of the galactocentric ra- 



dius R. The solid/dashed curves show the effect of unplug- 
ging/plugging the subhalo tidal disruption for both the Via 
Lactea II (red) and Aquarius (blue) setups. 



flux of a single object in a standard manner, with the 
factor 1 / y/N. A large value would express the fact that 
a small number of subhalos may have a large impact on 
the predictions, while a small value ensures that the pre- 
dictions are typified by contributions of a large number 
of objects. Of course, that depends on the volume rel- 
evant for the average, which is actually set by the CR 
propagation scale, beyond which a CR can hardly reach 
the Earth. Consequently, that volume is connected to the 
CR energy. For positrons, the high (low) energy range 
is characterized by short (long) propagation scales. The 
opposite is true for antiprotons. In that case, the aver- 
age number of subhalos expected in the corresponding 
small (large) volume is small (large) and is thus subject 
to large (small) fluctuations. Because the local density of 
the smooth DM component is rather low, a few objects 
can have a strong impact on the predictions, while the 
probability to find them in such a small volume is accord- 
ingly small. When the average volume, or equivalently 
the propagation scale, increases, the relative fluctuation 
is thus expected to decrease. Again, we stress that the 
energy dependence is reversed for antiprotons. Note fi- 
nally that the variance for positrons diverges like 1 /A 3 
in the limit of vanishing propagation scale and very local 
source term xs — > xq , so that we will have to put a cut- 
off to avoid numerical losses. This cut-off radius, namely 
the distance to the Earth below which we keep constant 
the subhalo population and disregard its fluctuation, can 
be reasonably set to 50 pc, without loss of predictivity. 



Indeed, increasing the cut-off radius to 100 pc does not 
affect our results significantly. 

The overall antimatter CR flux is the sum of the sub- 
halo component plus the smooth component, the latter 
being somewhat corrected for not carrying the whole DM 
mass anymore: 4> tot = (1 — /) 2 (/> sm +</>sub> where / is a cer- 
tain DM fraction that we discuss in more details below. 
Given the expressions of the smooth and clumpy contri- 
butions, the so-called boost factor is merely defined from 
their ratio as 



B(E) = (1 - ff 



sm (£) 



(27) 



This boost factor has to be understood as the expected 
average enhancement associated with the subhalo mod- 
eling. We emphasize that it does depend on energy, the 
ratio of the clumpy to smooth component having no rea- 
son to be constant with energy. Actually, as underlined 
above, the energy sets the propagation scale, and conse- 
quently the volume inside which the signal is integrated. 
The smooth component will dominate at large propaga- 
tion scale, when the smooth dense regions close to the 
Galactic center come into play, i.e. at low (high) energy 
for positrons (antiprotons). On the contrary, the contri- 
bution of subhalos can be significant at small propagation 
scales because the smooth contribution is more local in 
that case and hence much lower. There is a close parallel 
to make with the boost factor as computed for 7-rays, 
for which the relevant physical variable is not the energy, 
but the angle between the line-of-sight and the Galactic 
center direction ip: the boost is negligible at small angles 
because of the large contribution of the smooth central 
part of the Galaxy. Moreover, it is worth mentioning that 
the boost factor also depends on the injection spectrum 
for positrons. It is not the case at all for antiprotons, 
for which the injected spectrum cancels in the subhalo to 
smooth ratio (see end of Section VA for more details), 
though the associated boost factor still depends on the 
energy because of propagation. 

The fraction / appearing in the equation above is the 
fraction of dark matter in form of subhalos, that has to be 
subtracted to the smooth component: p sm — > (1 — /)/? S m 
(note that the definition of the smooth component in 
Eqs. [l] and [2] already accounts for this subtraction). As 
already stressed in [55] , this fraction depends on the way 
to normalize the problem. One can choose to keep the 
whole Galaxy mass constant when adding subhalos, in 
which case / — /m is the mass fraction in subhalos, 
or one can choose to keep the average local DM density 
constant, in which case / = /© is the local density frac- 
tion. These two normalization patterns are not equiv- 
alent, unless the spatial distribution of subhalos tracks 
the smooth profile. Therefore, if the Galaxy mass (local 
density, respectively) is kept constant, the average local 
density (Galaxy mass) will change after adding subhalos 
on top of a smooth component. For instance, we find 
that f M = 53.33%/17.50% and / = 2.96%/0.25% in 
our modelings of the Via Lactea II and Aquarius sim- 
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ulations, respectively, such that a normalization to the 
mass would decrease the average local density by a fac- 
tor (1 — /m) _1 ~ 2.0/1.2 compared to a smooth descrip- 
tion of the halo. Notice that although this is not that 
important for the computation of the boost factor, it 
could affect significantly the overall predicted fluxes at 
energies corresponding to small propagation scales, since 
in that case, the flux is cx (p© 2 ). For instance, if the 
mass normalization pattern is chosen for a Via Lactea II 
-like model, and if subhalos are subdominant in terms 
of flux (for large propagation scales), then one expects a 
decrease in the overall flux by a factor of ~ 2 2 = 4 with 
respect to a smooth description of the DM distribution. 
Keeping this issue in mind, we will take / = Jm in the 
following, to be consistent with the 7-ray analysis. 

The boost factor, as written above, stands for the av- 
erage enhancement of the DM signal in antimatter CRs 
due to the presence of subhalos in the Galaxy. It is as- 
sociated with a statistical fluctuation that is straightfor- 
wardly connected to that of the total subhalo flux: 

°e {E) " W) ■ (28) 

We see that even though fluctuations of the subhalo flux 
were found to be large compared to the subhalo flux itself, 
the boost factor would have a sizable variance only if 
those fluctuations are greater than the smooth flux. This 
mostly characterizes the large propagation scale regime, 
where the smooth contribution dominates the signal and 
completely overcomes the statistical variance expected 
from the subhalo flux. 

It turns out to be possible to derive an analytical ex- 
pression for the boost factor [55] , or equivalently for the 
total subhalo flux, in the vanishingly small propagation 
scale limit (i.e. at the very high/low energy tail of the 
positrons/antiprotons distributions). This asymptotic 
expression is very convenient not only to check numerical 
computations, but also because it usually corresponds to 
the maximal mean value of the boost factor — the boost 
can of course fluctuate around its mean value. This ana- 
lytical limit rests on the fact that at very short propaga- 
tion scale, the Green function G(xq x) 8(x — x@), 
so that Q —> S(x — af©)4g. We are therefore left with 
local quantities: 

B G = (l-fr + N^(C suh (R e )) Msub dP ^^(29) 

We emphasize that this expression is valid for any CR 
species and for any set of propagation parameters, in the 
regime of vanishingly small propagation scale. Never- 
theless, we also remind that such a regime is generally 
associated with large statistical fluctuations of the boost 
factor, because the average number of subhalos in such 
a small volume can be of the order of unity or less: the 
actual boost can be much larger if we sit on the top of 
a subhalo, or much lower if no bright object wanders in 
the neighborhood. 



C. Benchmark results and discussion 

Armed with the previous definitions and equations, we 
have computed the overall positron and antiproton fluxes 
and corresponding boost factors for the benchmark par- 
ticle physics scenarios defined earlier, and for the Via 
Lactea II and Aquarius cases. 



1. Boost factors 

In Fig. [lOj we show the results obtained for the boost 
factors and associated 1-a statistical bands, for both the 
positron and the antiproton signals. The top left panel 
represents the boost factor for positrons given different 
DM distributions (Via Lactea II- or Aquarius- like), as- 
suming the injection of a 100 GeV positron line at source. 

It is noteworthy that both DM subhalo models give 
similar and almost irrelevant enhancements, even less 
than unity when the smooth component dominates. In- 
deed, in this case, B « (1 — /) 2 , where / is the mass 
fraction in subhalos (see Eq. [27] and comments below) , 
which explains why the minimal value obtained for the 
Via Lactea //configuration is less than the one obtained 
with the Aquarius one (see the mass fractions in Table [I]). 
Nevertheless, the former still provides a slightly higher 
asymptotic boost factor than the latter, mostly because 
of the steeper mass function, which increases the relative 
contribution of the much more numerous and concen- 
trated lightest subhalos. This plot also shows that the 
effect of including tidal disruption in the model (dashed 
curves) is negligible except at low energy where it slightly 
reduces the boost factor. This behavior is due to the 
fact that tidal effects become relevant within a ~ 3 kpc 
from GC (or equivalently ~ 5 kpc from the Earth), corre- 
sponding to large propagation scales at which the smooth 
contribution dominates. 

In the middle panel of Fig. [lOj we show the effect of 
exploring the different benchmark particle physics mod- 
els introduced in Section IIIII The variation of the boost 
factor with energy looks less prominent for positron spec- 
tra than for a positron line, but this is mostly due to the 
fact that, for convenience, the x axis shows e = E/m x 
instead of the absolute energy. Had we found a larger 
asymptotic boost factor, the energy dependence would 
have been steeper. Notice that the variance is larger and 
more spread for differential injected spectra compared to 
the positron line. For the latter, indeed, positrons de- 
tected at lower and lower energies can originate from far- 
ther and farther regions, which smears the fluctuations; 
this is not the case for continuous injected spectra. 

Finally, in the right panel we show the boost factor 
obtained for antiprotons. This quantity does not depend 
on the benchmark model at all, since in this case the in- 
jected spectrum can be factorized out of the propagation 
description. The energy dependence comes only from 
propagation effects. The trend is the opposite of the 
positrons case. Indeed, small/large propagation scales 
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are found at low/high energy for antiprotons, contrarily 
to positrons. Note that the asymptotic values at vanish- 
ing propagation scales are exactly the same, as expected. 



2. Predictions of the antiproton and positron fluxes 



We have computed the expected antimatter flux in all 



which reads for standard quantities: 



the benchmark particle physics models of Section III 
The results are shown in Fig. |11[ where we have taken 
the solar modulation into account by applying a force 
field of 600 MV. In the left/right panel we show the 
positron/antiproton fluxes respectively, with the associ- 
ated 5cr fluctuation band due to the presence of subha- 
los. Theoretical predictions are compared with obser- 
vational data taken from [SB] [HS] HH] for positrons, and 
from [90, 91, 92, 93, 94J for antiprotons. We present the 
results for the Via Lactea //setup only since they are very 
similar to those obtained with the Aquarius setup. What 
will be quoted as smooth halo refers to a DM distribution 
without subhalos, while clumpy halo characterizes a DM 
distribution composed of a smooth component plus sub- 
halos. It is important to remind that the former is char- 
acterized by a scale density of p s = 0.31GeV/cm 3 , cor- 
responding to a local DM density p = 0.42GeV/cm 3 , 
while the latter has those parameters weighted by (1 — /), 
where / = 0.53 is the mass fraction in the form of subha- 
los (the smooth and clumpy halos have the same mass). 
Because the smooth DM density is reduced by a factor 
of ~ 2, we expect the flux predictions corresponding to 
the clumpy halo to be lowered by twice the same factor 
compared to those from a smooth halo at high/low en- 
ergy for antiprotons/positrons (large propagation scales). 
Likewise, at low/high energy for antiprotons/positrons, 
only the effect of subhalos can compensate for this loss 
in the clumpy halo configuration. From the calculation 
of the boost factor, shown in Fig. [10] it turns out that 
it is actually half the case, since the local boost is ~ 2. 
Therefore, quite interestingly, we end up with an aver- 
age prediction which is lower in the case of a clumpy 
halo than in the case of a smooth-only halo. This is the 
consequence of the normalization pattern adopted here, 
imposing the whole DM halo mass to be constant. That 
effect was already emphasized in |65j . 

Looking further at Fig. |11| only two benchmark mod- 
els appear for the antiproton flux, those with annihilation 
into bb and W + W~ . Indeed, the other final states con- 
sidered in this study do not cascade to antiprotons. 

Whereas our particle physics benchmark models are far 
from being observable in the antiproton spectrum, giv- 
ing primary fluxes well below the secondary background, 
positron signals could, at least in one case, exceed the sec- 
ondary background. Indeed, for the direct annihilation in 
e+e - , with m x — 100 GeV, we see that fluctuations due 
to subhalos are not needed to get a flux of the order of 
the secondary background. This can be understood eas- 
ily by deriving the general analytical expression of the 
flux for an injected positron line in the limit E — > m x , 
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Surprisingly enough, it is exactly the value of the pre- 
dicted background flux 



4> h e f(E 

at 100 



100 GeV) « 3 • 

lO-^cm^s^GeV^sr- 1 at 100 GeV in the median 
model of [58 . This formula is readily applied to the 
Via Lactea II setup by using p Q — 0.42GeV/cm 3 and 
multiplying by t he l ocal boost factor, which we can be 

100 GeV) « 



taken from Fig. 



We find (j) x + (my. 



5.5 • 10 cm s M^eV sr , which is larger than the 
background. We stress that this asymptotic flux predic- 
tion is only valid for E — > m x and falls thereby very 



quickly with m x like 



We will further comment 



on this when discussing the positron fraction in Sec- 
tion Em 

Apart from the positron line model, unfortunately, not 
even large fluctuations of the flux due to subhalos can 
lead to observable signals for the other benchmarks, as 
shown in Fig. [TT] This is in full agreement with the 
results obtained by [S3] from both a low resolution N- 
body simulation and analytical predictions. Of course, 
one could still invoke a single (or several) nearby mas- 
sive subhalo that could dominate the positron yield and 
provide a large enhancement, as early illustrated in |55j . 
even if the probability of such a configuration is van- 
ishingly small (it is actually accounted for in the calcu- 
lation of the boost factor fluctuation). Indeed, we do 
live in one realization of our Galaxy, which could corre- 
spond to this improbable configuration. In this case, if 
passing local kinematical constraints, a multimcssenger 
analysis using multiwavelength photons and antimatter 
cosmic rays could provide self-consistent predictions or 
constraints, which are quite easy to derive for one (or 
few) nearby object. Such a scenario may, however, al- 
ready be in tension with the EGRET data [53] , 

Finally, it should be noted that the predictions derived 
above rely on the median set of propagation parameters 
of [83 . The theoretical uncertainties on those parameters 
are still large, and the resulting uncertainty in terms of 
flux can easily reach one order of magnitude |83l I85j 



3. Comments on the positron fraction 

The excess in the positron fraction above a few GeV, 
recently made clearer with the release of PAMELA 
data PP, but previously hinted by the HEAT [55] and 
AMS data [97,, has triggered an impressive number of 
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FIG. 10: Left: boost factor on positrons for different (sub)halo models, Via Lactea II and Aquarius, assuming a 100 GeV 
positron line injected at source. Middle: boost for positrons for different injected spectra and for the Via Lactea //case. Right: 
boost for antiprotons, and for the same (sub)halo modes as in the left panel. All boosts are computed without solar modulation. 
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FIG. 11: Positron (left panel) and positron ratio (central panel) spectra for different injection spectra. Right: antiproton fluxes 
for bb and W + W~ spectra. All plots are done in the framework of the Via Lactea II simulation. 



studies, most of them dedicated to a possible DM inter- 
pretation (e.g. [98). Most of the predictions rely on the 
assumption that the DM annihilation rate is boosted, 
essentially from the non-relativistic Sommerfeld effect, 
with large branching ratios to leptons. Some others in- 
voke instead DM decay with a tuned lifetime. All of 
these assumptions are somewhat fine-tuned, and most of 
them do not treat the background consistently with the 
primary component. Here, we provide self-consistent pre- 
dictions for some benchmarks mostly motivated by parti- 
cle physics (except for one leptophilic model) and, more 
important, do not demand a good fit to the PAMELA 
data. Instead, we aim at testing the potential imprints 
and the detectability of such scenarios in the antimatter 
spectrum. 

In the central panel of Fig. [TT] we plot the results for 
the positron flux in terms of the corresponding positron 
fraction, i.e. 4> e +l(4>e+ + 4>e-)- The data appearing on 
the plot have been taken from [H [9GJ |97l |99] . The actual 
denominator of the positron fraction can be obtained by 
fitting the Fermi data on the sum of cosmic electrons plus 
positrons [TH] above 20 GeV, avoiding thereby to invoke 



any model of primary or/ and secondary component. At 
lower energy, we have constrained the electron spectrum 
from a fit on the AMS data [50] , and for positrons, we 
have taken the sum of our primaries plus secondaries, 
assuming that there are no other sources of positrons. 
That assumption is conservative in the sense that below 
10 GeV, the secondaries dominate over all our primaries. 
We have linked the two domains by interpolating over a 
range of a few GeV. 

As already discussed in Section |V C 2| only the WIMP 
model annihilating in e + e~, with m x — 100 GeV, pro- 
vides a sizable contribution to the positron fraction, over- 
topping the secondary background above 10 GeV. The 
other benchmark models contribute at most at the per- 
cent level, even when considering 5a subhalo fluctua- 
tions. Of course, our leptophilic model could afford for 
a large part of the PAMELA data. Anyway, we remind 
the reader that viable astrophysical explanations exist, 
like the contribution of local pulsars [75], and it might 
be difficult to distinguish between those different solu- 
tions to the positron excess, given the limited sensitivities 
and energy resolutions of the current experiments (even 
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a x 2 analysis is hardly relevant due to the amount of 
theoretical uncertainties). Our leptophilic model would 
give a sharp cut-off in the positron fraction as well as 
in the positron flux above 100 GeV, which unfortunately 
would appear smeared because of energy resolution ef- 
fects. Positron data at higher energy would really help 
to clarify this issue of the possible contribution of DM 
annihilation to the positron spectrum. Indeed, if the ex- 
cess is still prominent above say 200 GeV, a leptophilic 
model could hardly fuel the dominant contribution with 
conventional parameters, due to the rn~ 4 scaling of the 

-3.5 



primary flux, when E - 
scaling of the secondary background (see Eq. 30 ) 



compared to the E 



VI. CONCLUSIONS 

We have studied the 7-ray and antimatter fluxes aris- 
ing in two astrophysical setups built on the recent Aquar- 
ius and Via Lactea II high-resolution N-body simula- 
tions. Our aim was to quantify the differences between 
the two simulations in terms of prospects for detection, 
and to assess the impact of extrapolating the mass func- 
tion and concentration of subhalos down to their min- 
imum mass, that for common DM candidates can be 
10 order of magnitude smaller than the mass resolution 
of simulations. We stress that the two simulations are 
in remarkable agreement with each other, and that the 
biggest differences among in the prospects for detection 
arise from the different extrapolations of physical quan- 
tities suggested by different groups of authors (as in the 
case of the mass function of subhalos), and from the ap- 
plication of the results of the simulations to the specific 
MW halo. We present our results for the Aquarius and 
Via Lactea II setups, meaning that they correspond to 
models that we have built on simulations, and not to the 
simulations themselves. 

We have studied the different contributions to the 7- 
ray flux arising from the smooth DM halo of the MW, 
resolved subhalos (that we have generated with a Monte 
Carlo procedure), unresolved subhalos and extragalactic 
halos and subhalos. The smooth component dominates 
the annihilation flux in the inner regions of the galaxy. 
In the Via Lactea II setup, the resolved and unresolved 
subhalos dominate the annihilation flux at angles larger 
than ~ 20 degrees from the GC, all the way to the anti- 
center, where the extragalactic flux becomes comparable, 
though never dominant. In the Aquarius setup, the sub- 
structures component is suppressed, and never exceeds 
the smooth component. The extragalactic flux becomes 
instead dominant at angles larger than ~ 60 degrees from 
the GC. 

We have provided full sky maps that can be used as 
templates for DM searches with current experiments such 
as Fermi. If the search is concentrated towards the GC, 
there is little difference between the two simulations, 
in the sense that the profile of the annihilation flux is 
very similar in the two setups, while the normalization 



must be kept free given the uncertainties on the particle 
physics parameters. The optimal strategy to search for 
an annihilation signal is to take the sky maps provided 
by Fermi, that will constitute the 'background', and take 
our DM templates as signal. One can then easily esti- 
mate the size and shape of the region around the GC 
that maximize the S/N ratio. 

For a fixed particle physics model, the annihilation flux 
from the central regions of the Galaxy in the Via Lactea 
II setup is smaller by a factor O (10) with respect to 
Aquarius. This is due to the i) the larger fraction of mass 
in substructures in the Via Lactea II setup, that reduces 
the annihilation flux from the smooth component by a 
factor (1 - fc[ q ) 2 /(l - fYi L2 ) 2 ~ 3; 11) the smaller local 
density, i.e. DM density in the solar neighborhood, in the 
Via Lactea 77 set up, that reduces the the annihilation flux 
from the smooth component by a factor (/? Q 9 / p^ ) L2 ) 2 = 
(0.48/0.2) 2 = 5.8 and Hi) the smaller total mass. If one 
rescales the Via Lactea II and Aquarius setup to match 
the most recent determinations of the local density, and if 
the same subhalo mass fraction f c i = 0.18 is adopted for 
both simulations, then the annihilation maps look almost 
identical. 

Should the search for the diffuse emission from the 
GC fail because of the complicated astrophysical back- 
grounds in what is probably the most crowded region of 
the sky, the possibility remains to search for unidentified 
gamma-ray sources, that would appear as non-variable 
bright spots with no astrophysical counterpart, possibly 
correlated with dwarf galaxies, and with identical spec- 
tra. The number of detectable sources in both simulation 
setups is very similar, and for an optimistic DM scenario 
is between 1 and 10 for the Fermi-LAT in 5 years of op- 
eration. 

Finally, we have calculated the antimatter fluxes in 
both simulation setups, and we found that the boost fac- 
tor often invoked to provide a viable DM intepretation 
of the cosmic leptons puzzle, are completely unrealistic. 
The only annihilation channel that provides a sizeable 
enhancement of the positron ratio is direct annihilation 
to e + e~, which provides a flux higher than the secondary 
background for the set of propagation parameters used 
here, even without the help of any subhalo contribution. 
Although this model seems an interesting possibility, we 
want to stress that (i) it has been tuned to provide an 
exception case to the usual need of large boost factor 
to interprete the PAMELA data (see Eq. 30 and com- 



ments below) without any particle physics motivation, 
(ii) slightly increasing the DM particle mass above 100 
GeV would completely erase such a peak with respect to 
the background because of the dependency of the peak 
amplitude in rn x ~ A and the steep decrease of the flux at 
lower energies, (iii) we did not include the contributions 
of other astrophysical primary sources and (iv) the fit to 
the PAMELA data is still rather poor. 
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